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ABSTRACT 


A general approach is developed for predicting the power output of a concentrator 
enhanced photovoltaic space array. A ray trace routine determines the concentrated 
intensity arriving at each solar cell. An iterative calculation determines the 
cell's operating temperature since cell temperature and cell efficiency are functions 
of one another. The end result of the iterative calculation is that the individual 
cell's power output is determined as a function of temperature and intensity. Circuit 
output is predicted by combining the individual cell outputs using the single diode 
model of a solar cell. Concentrated array characteristics such as uniformity of 
intensity and operating temperature at various points across the array are examined 
using computer modeling techniques. An illustrative example is given showing how 
the output of an array can be enhanced using solar concentration techniques. 
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INTRODUCTION 


The incentive for using a concentrated array is to improve power to mass 
performance and cost effectiveness. Many different schemes for providing an array 
with concentrated illumination have been proposed. Systems using lenses, mirrors, 
or combinations of mirrors and lenses have been suggested. Finding the optimal 
concentration method for any particular situation (near earth or deep space for 
example) represents a challenging engineering problem. 

The only certain way of determining the performance of a particular concentratioi 
method is to actually build a scale model using proposed materials and solar cells. 
However, the high cost associated with building and testing such scale models makes 
it more reasonable to first model the power performance of various concentration 
concepts using computer techniques. Such modeling can dramatically narrow the 
choice of potential concentration schemes by discarding the least favorable candidates 
While the modeling approach is general, the modeling results are highly specific, 
based upon the particular selection of solar cell type, circuit size, circuit 
orientation and concentrator geometry. 

The concentrated array selected for illustration is the compound parabolic 
concentrator (CPC) in combination with a 2-ohm-cm, 2x2 cm, 50-micron thin silicon 
solar cell array. The CPC consists of truncated parabolic specular mirrors extending 
along two sides of the rectangular array. (Reference 1 and 2). The geometric 
intensity concentration ratio (Cg) can be adjusted by changing the angle of the 
mirrors to the array (Figure 1). The array selected consists of a single circuit 
of four cells in parallel by 20 cells in series. Rays entering the concentrator 
either strike the array directly or are redirected by the mirrors onto the array. 

The first step in the modeling process is to trace a sufficient number of evenly 
spaced incident rays to determine the intensity profile across the array. Each 
successive modeling step, as it applies to this illustrative example, is examined 
in the subsequent sections. 
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CONCENTRATED INTENSITY DETERMINATION (RAY TRACE ROUTINE) 

One thousand evenly spaced rays enter the CPC's entrance aperture and are traced 
until they encounter the array. With the CPC mirrors aligned in the maximum con- 
centration configuration (Figure 2), all rays experience at most a single reflection 
before striking the array (assuming perfectly specular mirrors). As the concentration 
ratio (Cg) is decreased by tilting the tips of the mirrors inwards through an angle 
theta, multiple reflections are possible. However, all rays still reach the array. 

If mirror imperfections or Intentionally stippled mirrors are considered (stippled 
mirrors are useful in reducing intensity nonuniformities), rays can be reflected 
back out the concentrator's entrance aperture. Such rays do not contribute to the 
generation of electrical power. The illustrative example deals exclusively with 
the specular mirror case. 

The power per ray is determined from equation (1) 

Total incident power entering concentraotr's aperture 
Power of single ray = arbitrary length along the axis of symmetry 

Total number of rays 


_ lAC 
T0OO 

where I = incident intensity at concentrator's aperture (W/m2) 

A = width dimension of concentrator aperture (see Figure 1) (m) 

C = arbitrary length along axis of symmetry (m) 

The power per ray after a reflection is given by equation (2): 

Power (postreflection) = Power (prereflection) • (1-R) (2) 

where R is the mirror reflectivity, (For a discussion of mirror reflectivity 
see Appendix A). 
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ENTRANCE APERTURE A 



ARRAYS LOCATED AT 
CONCENTRATOR' S EXIT APERTURE 


Figure 2. Typical Rays Traced through a CPC at the Haximum 
Concentration Configuration 
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The CPC exit aperture (colocated with the array) Is divided Into 50 bins for 
the purpose of counting rays. The Incident power on each bln Is the sum of the 
power of the 1nc1d^»nt rays. The Intensity for each bln 1s the power per bln divided 
by the bln area [bln area = (B-C)/50]. B * width dimension of array (see Figure 1) (m). 
Note that the bln length (C) cancels Itself In bln Intensity and hence Is an arbitrary 
value. 


TABLE 1 

Input Parameters for a Variable CPC Enhanced Array 


6 

(degrees) 

Cg 

Ml rror 
height 

(cm) 

"A" 

width 

(cm) 

ml rror 

reflectivity * 

0 

4.63 

35.76 

37.73 

.865 

1.9 

4.34 

36.23 

35.34 

.860 

3.9 

4.02 

36.68 

32.80 

.860 

6.0 

3.69 

37.11 

30.90 

.855 

8.2 

3.34 

37.50 

27.22 

.850 

10.3 

3.00 

37.83 

24.46 

.845 

13.5 

2.48 

38.22 

20.22 

.850 

17.9 

1.76 

38.58 

14.31 

.900 

* see Appendix A 
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The first and the fiftieth bins are disregarded since cells are not positioned 
immediately next to the base of the mirrors. Rays entering these uncelled bins do 
not generate electrical power. The remaining 48 bins are assumed to overlay the 
4 solar cells (see Figures 1 and 2). The intensities of the 12 bins overlaying 
each cell are averaged together to determine the average intensity reaching each solar 
cell. Table 1 shows how the geometric concentration ratio (Cg = A/B) and other 
parameters vary as the mirrors are tilted inwards through the theta angle measured 
from the maximum concentration configuration. 

If the concentrator has a twist along the axis of symmetry, the intensity 
profile calculation must be repeated at various stations along the length of the 
concentrator. However, since this is not fundamentally more difficult, the no- 
twist case is illustrated. 

The model also has the capability to introduce distortions into the mirrored 
surfaces. These can be either random in nature, simulating imperfections in specular 
surfaces, or periodic, simulating intentional mirror patterning, useful in reducing 
intensity nonuniformities. The illustrative example deals exclusively with the 
zero distortion, specular mirror case. 

SOLAR CELL OPERATING TEMPERATURE DETERMINATION 

Assuming that the solar cell is in good thermal contact with a substrate of 
known absorptivity and er.ilssivity, the temperature of a solar cell operating in 
space can be calculated from Stefan's Law: 


t 4 , («-n(T))S 


(3) 
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where 


T = cell operating temperature (K) 

a “ absorptivity (dimensionless) 

n(T) ■ solar cell efficiency (itself a function of 
temperature and intensity) (dimensionless) 

S = concentrated solar intensity (W/m ) 

= emissitivity of front surface (dimensionless) 

= emissitivity of back surface (dimensionless) 

-8 2 4 

a * Stefan-Boltzmann Constant = 5.6696 • 10 W/m K 

Since n» cell efficiency, is a function of temperature, it is necessary 
to assume an initial value of n and solve equation (1) by iteration. 

This iterative process converges rapidly but a computer based calculation 
avoids tiresome repetitive calculations and interpolations. Cell 
efficiency is interpolated from tabulated data as explained in the 
following section. 

INDIVIDUAL CELL OUTPUT DETERMINATION 

The individual cells short circuit current (Isc), maximum 
power current (Imp), open circuit voltage (Voc) and maximum power 
voltage (Vmp) are determined from tables of experimentally measured 
cell data (Reference 3). The tabular data set for each current and 
voltage is given between temperature and concentrated intensity limits 
of -160°C to 140°C and 5 mW/cm to 250 mW/cm, respectively. Data points 
are given at 20®C increments in temperature and at four intermediate 
intensities. Given an operating cell temperature and intensity, a 
linearly interpolated current or voltage is determined. Figures 3 and 
4 show the respective variation of Voc and Vmp with intensity. While 
logarithmic interpolation is more appropriate than linear interpolation 
for voltages, the difference between the two interpolation methods 
typically amounts to less than 1 mV. This is well below the measured 
data standard deviation for either Voc or Vmp. Hence, the linear inter- 
polation method was adopted for both current and voltage because it 
simplified the interpolation calculations with no meaningful loss in 
accuracy. 

A limitation of the existing data is that the upper intensity 
limit of 250 mW/cm is exceeded for the illustrative example. The 
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-olarex 50 micron 2a -cm Voc (mV) 








Figure 4. 2 x 2 cm Solarex Micron Vmp 







existing data has been graphically extended to encompass these higher Intensity 
regions (dashed lines of Figures 3 and 4). Experimental confirmation of these 
graphical extensions would be highly desirable. 

PARALLELED CELL STRING OUTPUT DETERMINATION 

Paralleled cells In a compound parabolic concentrator enhanced array experience 
widely different temperature/ Intensity conditions. Localized intensities of over 
1250 mW/cm (over 9 suns) at an Incident intc^nslty of 67.65 mW/cm (half sun) are 
predicted. Under these conditions, localized temperature hot spots reach 275°C. 

The resulting cell's current and voltage outputs also vary widely. Similar, though 
less extreme variations, can be caused by Inherent cell-to-cell differences. The 
resulting mismatches In cell powers reduces the total available power below that of 
the simple sum of Individual cell powers. The single diode model of a scicn' cell 
Is used to predict the output of the mismatched paralleled cells. This technique 
essentially generates an I-V curve for each paralleled cell. The single diode 
model parameters for each cell are adjusted so that they agree with the I sc. Imp, 
Voc, Vmp, and series resistance (R^) values of each cell In the parallel string. 

Cnee the single diode model parameters are made to match the observed values of 
each cell, then the cell's entire I-V curve Is uniquely specified. The Individual 
\ cell's I-V curves may then be algebraically added to give the total paralleled 

string I-V performance curve. From this curve, the string's outputed power Is 
easily determined (Reference 4). 

The equations necessary to determine the Individual cell I-V curves will now 
be derived from the basic single diode equivalent circuit for a solar cell, which 
follows: 



WHERE IL IS A CONSTANT 
CURRENT GENERATOR 
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Let: 


r “load current 
c 

= load voltage 

I = reverse saturation current 
0 

I|_ = light generated current (a constant current generator) 

6 h “ shunt conductance (an initial estimate provides for 
^ faster convergence but is not required) 

= diode current 

c - MI 

0 ■ q 

A = an assumed constant 
k = Boltzmann constant 
T = temperature (K) 
q = electron charge 
Rg = series resistance 
I sc = short circuit current 
Imp = maximum power current 
Voc = open circuit voltage 
Vmp = maximum power voltage 
Pmax = Imp • Vmp = maximum power 

From the basic diode equation: 

II lo [exp -1] + (V(, + I(. Rg) Gsh + (4) 


Evaluating (2) at the short and open circuit conditions one obtains respectively: 
IL = lo Cexp -1] + Isc (1 + Rs Gsh) (5) 


and 

IL = lo [exp (|^^) -1] + Voc Gsh (6) 


il 


At the maximum power point 


IL = lo [ exp (!!!!E-+J!!!E_Si) .]] + Vmp Gsh + Imp (1 * Rs Gsh) 
dP 

also since ^ = 0 at the maximum power point; 


diffe'^entlatlng (2) with respect to V^. and using (6), one obtains: 


M ^sh^feexp 

Vmp 1 ^ p, Gsh ^ ye. e'p ( Vmp . ^^ImO Ii' 

Using (3), (4), (5) and (7) to eliminate Il> Ushe and Iq* one obtains; 


(7) 


( 8 ) 


(9) 


En = 


(Vmp-Imp Rs) (Voc-Isc Rs) (Isc Vmp-Isc Voc Imp Voc) 

^ Vmp (Voc-Isc Rs)(2Imp-Isc)-(Voc-Isc Rs)(Imp Voc-Isc Vmp)- ^ Imp(2Vmp-Voc)(Voc-Isc ^ 

A(n "in 


no) 


where 

Xs = exp Xo - exp (|^); X„ = exp ( ) 


Equation (10) Is solvable by iteration. 

Once Eq is known then Iq and Gsh are obtained from; 

T _ Isc Vmp-Isc Voc Imp Voc 

° " (Xo-Xs) (Vmp + Imp Rs)-(Xm-Xs) Voc-Isc Rs (Xo-Xm) 

and 


*>sh = llr -■ 


Voc - Isc Rg 

IL Is then obtained from (4); Ic as a function of V^. is then obtained (I.e. the I-V 
curve Is obtained). In this manner, the Ij. of each paralleled cell can be calculated for 
every Vc. Summing the Individual cell Ic values for given values of Vc yields the total 
paralleled string I-V curve. 
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ARRAY OUTPUT DETERMINATION 


As currently configured* the model assumes that the array consists of a single 
circuit. The model can be generalized to an array consisting of an arbitrary set 
of n circuits. As configured, the I-V curve for the circuit Is the algebraic sum 
of the paralleled string curves. 

Table 2 summarizes the Appendix B results for the variable 
concentration CPC enhanced thin cell array with specular mirrors. The 
theta angle Is the angle the mirrors have been tilted Inwards from 
the maximum geometric Intensity design point. The figure of merit W 
shown In the last column of Table 2 Indicates the power gain of the 
concentrator enhanced thin cell array as referenced to an Identical 
unconcentrated thin cell array (identical thermal properties) operating 
under identical incident Intensity conditions (see Fig. 5). 


TABLE 2 

Summary of Appendix B Results for the 
Variable CPC Thin Cell Array with Specular Mirrors* 


mirror angle 

geometric intensity 

predicted 80 cell array module output 

"W" 

0 (degrees)** 

concentration 

ratio 

voltage(V) 

amperage (A) 

power(W) _ 





Vmp 

Imp 

PMP 


0.0 

4.629 


2.200 

.804 

1.769 

.66 

1.9 

4.336 


2.600 

.808 

2.100 

.79 

3.9 

4.024 


4.000 

.792 

3.166 

1.19 

6.0 

3.692 


5.600 

.752 

4.209 

1.58 

8.2 

3.340 


6.400 

.711 

4.551 

1.71 

10.3 

3.002 


6.000 

.615 

3.692 

1.38 

13.5 

2.480 


6.200 

.506 

3.138 

1.18 

17.9 

1.756 


7.800 

.397 

3.094 

1.16 

*Assumpt1ons as follows: 

80 cell 

module, 2x2 

cm cells, 4 

cells In parallel 



spread across width of the CPC; incident unconcentrated 
intensity 67.65 mW/cm^; array thermal absorptivity 
= 0.84, array front hemispherical emissitivity = 0.83, 
array back hemispherical emissitivity = 0.76; cell 
series resistance = .2 ohms. 

**Mirror angle referenced to the maximum geometric intensity concentration 
design point. 

Figure 5 is a smoothed graphical representation of Table 2. 
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THETA 


Figure 5. The CPC Power Gain (W) vs Mirror Angle Tilt 
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SUMMARY OF RESULTS FOR THE ILLUSTRATIVE EXAMPLE 

The CPC enhanced thin cell array with specular mirrors does not 
generate Its peak power at the maximum intensity design point. The 
dominant cause is that the specular CPC is plagued by nonuniform 
intensity patterns which produce localized intensities as high as 
1250 mW/cm^ (over 9 suns) at an incident intensity of 67.65 mW/cm^. 
These nonuniformities result in localized hot spots which reach 
275“C. Maximum power was generated with the mirrors angled inwards 
8.2 degrees from the maximum intensity design point. While this 
concentrator design is almost certainly not the optimal design for 
this operating environment, the power gain as referenced to 
an unconcentrated array, operating under similar incident intensity 
conditions, is a factor of 1.7. Whether such a power gain is 
sufficient to pay back the added cost, mass, and complexity of the 
concentrator system remains an open question. 

CONCLUSION 

The general computer modeling approach developed in this memo 
provides a method for predicting the power output, uniformity of 
intensity, and operating temperature of a concentrator enhanced solar 
array. This modeling technique can be used to evaluate the relative 
merits of various concentrator concepts. Based upon the results of 
additional modeling, candidate concepts can be proposed for scale model 

testing. 
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APPENDIX A 


REFLECTIVITY MODEL FOR MIRRORED SURFACES 


APPENDIX A: Reflectivity Model for Mirrored Surfaces 


The current model provides mirror reflectivity as a function of 
angle of incidence. 


TABLE A-1. Mirror Reflectivity vs Angle of Incidence 


angle of ^icidence 
(degrees) 

mirror reflectivity 

0 (normal to mirror surface) 

.90 

5 

.90 

10 

.90 

15 

.90 

20 

.90 

25 

.90 

30 

.89 

35 

.89 

40 

.89 

45 

.89 

50 

.89 

55 

.88 

60 

.88 

65 

.87 

70 

.86 

75 

.85 

80 

.84 

85 

.87 

90 

1.00 


The'^ngle of incidence" here refers to the angle an incident ray makes 
to the plane connecting the extreme inboard and outboard edges of the 
mirror. While a single value of reflectivity does not allow for the 
changing angle of incidence along the curvature of the mirror, it does 
represent an "average" value of reflectivity. 
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Reflectivity values are based upon Information supplied by 
Dr. B. Zeldin. Dr. Zeldin supplied a matrix of reflectivity values 
for an aluminized mirror at a specific wavelength and a specific 
angle of Incidence (see Appendix A: Table ^~ 2 ). A single value of 
reflectivity for each angle of Incidence Is obtained by weighting 
the reflectivity at a particular wavelength by the amount of energy 
present In the solar spectrum at that particular wavelength. Table 
A-3 of Appendix A contains the weighting factors for each wavelength. 
The final result Is the angle dependent reflectivity shown In 
Appendix A, Table A-1. 
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ty of Mirrored Surfaces as i Function of Angle of Incidence and Waveguide 
ed by Dr. B. Zeldin in private communication) 


ANGLE OF INCIDENCE (DEGREES) 



5 


15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 

80 

85 

90 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.90 

.90 

.90 

.90 

.89 

.89 

.89 

.90 

.94 

1.00 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.90 

.90 

.90 

.90 

.90 

.89 

.89 

.88 

.88 

.89 

.92 

1.00 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.89 

.89 

.88 

.87 

.87 

.88 

.91 

1.00 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.89 

.89 

.88 

.87 

.87 

.86 

.90 

1.00 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.89 

.89 

.88 

.87 

.86 

.85 

.88 

1.00 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.90 

.89 

.89 

.89 

.88 

.87 

.86 

.85 

.83 

.85 

1.00 

.89 

.89 

.89 

.89 

.89 

.89 

.89 

.89 

.89 

.88 

.88 

.88 

.87 

.86 

.85 

.83 

.81 

.83 

1.00 

.88 

.88 

.88 

.88 

.88 

.88 

.88 

.88 

.87 

.87 

.87 

.86 

.86 

.85 

.83 

.81 

.80 

.81 

1.00 

.87 

.87 

.87 

.87 

.87 

.87 

.87 

.87 

.87 

.87 

.86 

.86 

.85 

.84 

.82 

.80 

.79 

.80 

1.00 

.86 

.86 

.86 

.86 

.86 

.86 

.86 

.86 

.86 

.86 

.85 

.85 

.84 

.83 

.81 

.79 

.77 

.79 

1.00 

00 

.88 

.88 

.88 

.88 

.88 

.88 

.88 

.88 

.88 

.88 

.87 

.86 

.85 

.84 

.82 

.80 

.81 

1.00 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.91 

.90 

.90 

.89 

• 

00 

.86 

.83 

.83 

1.00 




TABLE A-3. Weighting Factor for Mirror Reflectivity vs. Wavelength 


Wavelength (A) 

Weighting Factor 

4000 

.71 

4500 

1.00 

5000 

.97 

j500 

.86 

6000 

.83 

6500 

.76 

7000 

.68 

7500 

.62 

8000 

.55 

8500 

.50 

9000 

.44 

9500 

.41 
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appendix b 


RESULTS OF THE VARIABLE CPC WITH SPECULAR MIRRORS 
IN COMBINATION WITH THE 50-MICRON SOLAREX CELL 
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appendix B contains a graphical sunwiary of the intensity, temperature and 
I-V curve results for each mirror tilt angle. Each of the temperature and intensity 
plots gives a graphical representation of temperature and intensity variation 
across the array width. The I-V graphs show three I-V curves per graph. The 
lower curves are individual cell curves while the upper curve is the algebraic 
sum of the individual cell curves. Only two of the four individual cell curves 
need be considered since the other two possess identical characteristics because 
of concentrator symmetry. Some of the I-V curves have been truncated by computer 
limitations, however, the curves still demonstrate the circuit effect of intensity/ 
temperature induced cell mismatches. 
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PNGLE 0F INCIDENT RHDIHTION IS 0.00 DEGREES RELATIVE TO OPTIC AXIS 



Figure B-1. Intensity vs Bin No. for Theta = 0.0 deg Configuration 


ANGLE 0F INCIDENT RADIRTI0N IS 0.00 DEGREES RELATIVE TO 0PTIC AXIS 



Figure B-2. Temperature vs Bln No. for Theta = 0.0 deg Configuration 





Figure B-3. Array I-V Curve for Theta « 0.0 deg 


fiNGLE 0F INCIDENT RROIRTIDN IS 0.00 DEGREES RELfiTIVE TO OPTIC RXIS 



Figure B-4. Intensity vs Bin No. for Theta * 1 .9 deg Configuration 
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ANCLE 0F INCIDENT RPDIRTION IS 0.00 DEGREES RELATIVE TO OPTIC AXIS 



Figure B-7. Intensity vs Bln No. for Theta = 3.9 deg Configuration 







Figure B-9. Array I-V Curve for Theta = 3.9 deg 



RN&LE OF INCIDENT RRDlfiTlON IS 0.00 DEGREES RELRTIVE TO OPTIC RXIS 
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Figure B-10. Intensity vs Bin No. for Theta = 6.0 deg Configuration 


flNGLF, nF INCH3ENT KRDIRTIGN IS 0-00 DEGREES RELATIVE TO OPTIC AXIS 



Figure B-71. Ten^rerature vs Bin No. for Theta = 6.0 deg Configuration 




Figure B-12. Arr^y I-V Curve for Theita 6.0 deg 
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Figure B-13. Intensity vs Bln No. for Theta » 8.2 deg Configuration 




RNGLE 0F INCIDENT RfiOIflTION IS 0-00 DEGREES RELATIVE TO ePTIC ftXI 



Figure B-14. Temperature vs Bin No. for Theta » 8,2 deg Configuration 





RNGLE 8F INCIDENT RHDiflTION IS 0-00 DEGREES RELRTIVE TO OPTIC AXIS 



Figure B-15. Array I-V Curve for Theta = 8.2 deg 


PNGLE OF INCIDENT RRDIRTION IS Q.QO DEGREES RELRTIVE TO OPTIC RXIS 



Figure B-16. Intensity vs Bin No. for Theta = 10.3 deg Configuration 



ANGLE 0F INCIDENT RADIATION IS 0.00 DEGREES RELATIVE TO OPTIC AXIS 



Figure B-17, Temperature vs Bin No. for Theta = 10.3 deg Configuration 





Figure B-18. Arruy I-V Curve for Theta » 10.3 deg 



PNGLE OF INCIDENT RHDIRTION IS 0-00 DEGREES RELPTIVE TO OPTIC PXIS 





Figure B-19. Intensity vs Bin No. for Theta = 13.5 deg Configuration 
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RNGLE OF INCIDENT RRDIRTION IS 0.00 DEGREES RELRTIVE TO OPTIC AXIS 
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Figure B-22. Intensity vs Bin No. for Theta * 17.9 deg Configuration 
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